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Abstract 



The spectral stabihty of hquid metal differentially rotating in trans- 
verse magnetic field is studied numerically by solving the eigenvalue 

^ • problem with rigid- wall boundary conditions. The equilibrium veloc- 

ity profile used in calculations corresponds to the electrically driven 
flow in circular channel with the rotation law 0,{r) oc 1/r^. This type 
of flow profile is planned to be used in new experimental device to test 
the magnetorotational instability (MRI) in laboratory. Our analysis 

QQ I includes calculations of the eigen-frequency spectra for both axisym- 

metric and non-axisymmetric modes. It is found that for chosen device 
parameters the flow is always spectrally unstable due to MRI with the 

^^ ' fastest growth rate corresponding to the axisymmetric mode. 

o' 
Cu\ 1 Introduction 

6, 

^ , Magnetorotational instability (MRI) - the instability of conducting fluid, 

j^ I rotating in magnetic field - is one of the most important processes in astro- 

physics. MRI is believed to be responsible for angular momentum transport 
in such objects as accretion disks formed around massive black holes and 
active galactic nuclei. This instability was originally discovered in 1959 by 
Velikhov [L^, who, though, did not consider it in an astrophysical context. 
It was suggested only in 1991 |2| that a weak magnetic field can linearly 
destabilize the otherwise stable accretion disk thus directly leading to disk 
turbulence and orbital angular momentum transport. MRI becomes opera- 
tive when the angular velocity decreases outwards, dQ/dr < 0, which holds 
for Keplerian flow in accretion disks (f2 oc r~^/^). 

Because of its importance in understanding of accretion disk physics, 
a number of experiments have been proposed recently |S|-|H] to study this 
instability in the laboratory. In most of them a Couette flow of liquid metal 
between two rotating cylinders placed in external magnetic field is supposed 
to be used However this simple design has its disadvantages. The main one 



is that the equihbrium flow is controlled almost entirely by the endcaps of 
the vessel - so-called Ekman effect Pj. In this case the velocity profile is 
significantly different from Couette flow, so the conditions for MRI may not 
be met. 

Another way to rotate liquid metal is to apply radial electric current be- 
tween cylinders placed in vertical magnetic field, that produces a Lorentz 
force driving metal. It is worthwhile noting that in this case the effect of the 
Ekman layers at the top and bottom insulating caps is considerably reduced. 
Actually one has to take into account only Hartmann layers which are very 
small for large Hartmann numbers and do not affect the main body of the 
flow. This idea was proposed in Russian Research Center "Kurchatov Insti- 
tute" by E. P. Velikhov (2004, private communication) and the experimental 
device is being designed there. A similar setup using the plasma instead of 
liquid metal has been developed in Los Alamos fij. 

In this paper we present a linear stability analysis of the electrically driven 
liquid metal flow in the circular channel for the configuration relevant to 
Kurchatov experiment. It has been shown in jTUj that in the stationary 
state the angular velocity profile of such flow scales with radius as fi oc 
1/r^ virtually throughout the entire channel cross section if the Hartmann 
number is large enough (Ha^ 1). This dependence is a particular case of 
the more general ones, Q = a + h/r"^ (Couette flow) and VL oc r~^ (with 
q = 3/2 corresponding to Keplerian flow) whose stability has been studied 
analytically and numerically in numerous papers. Therefore some results 
obtained in those papers are relevant to our study. Let us briefly review 
them. 

Originally MRI was discovered as a powerful local axisymmetric (with 
azimuthal wave-number m = 0) instability ^1 |2]. Consideration of non- 
axisymmetric modes {m ^ 0) in local approach (for example, JT]) has shown 
that these modes become stabilized with increasing values of m. 

Importance of global modes and boundary conditions for MRI problem 
have been first analyzed in ^j. In ^^ authors considered in details non- 
axisymmetric global modes with free boundary conditions which are not ap- 
propriate for the system with rigid walls. 

A great deal of work in global stability study of magnetized Couette flow 
(including our case with a = 0) in the frame of non-ideal magnetohydro- 
dynamics (MHD) has been done by Riidiger et al. (see, for example, |3], 
jH] and Jl]). They found marginal stability for different parameters of the 
system, but they did not analyze the growth rate of MRI and its dependence 
on azimuthal wave-number m. 

The goal of this paper is the global spectral stability analysis of the in- 
compressible MHD flow with angular velocity profile Q oc 1/r^. The eigen- 



value problem for this configuration is set up in Section 2. The main body 
of results were obtained numerically and are presented in Sections 3, 4, 5. 
In Section 3 axisymmetric modes {m = 0) are considered and their clas- 
sification according to radial wave-number is made. Section 4 deals with 
non-axisymmetric modes (m 7^ 0), the dependence of MRI growth rates on 
azimuthal m and axial k wave-numbers is presented there. In Section 5 we 
study marginal stability of the system under consideration. In Section 6 we 
discuss the possibility of MRI in electrically driven liquid metal fiow. 

2 Statement of the Problem 

We consider the following model of the experimental device proposed in Kur- 
chatov Institute (Fig. [T]). A liquid metal (Sodium) occupying a circular 
channel is driven by applying uniform vertical magnetic field Bq and ra- 
dial electric current Jq. The stationary state of this configuration has been 
studied in details in [TU]. It has been shown that in the cylindric system of 
coordinates (r, ip, z) the fiuid rotates around the 2;-axis with the angular 
velocity, which has the form 

f](r) = ^ (1) 

almost entirely in the channel if the Hartmann number is large (Ha^ 1). 

Despite the fact that the value of stationary velocity is determined by the 
dissipative effects (conductivity a and viscosity z/ of the fluid) in our study we 
neglect all of them and assume that the liquid metal is a perfectly conducting, 
ideal, incompressible fluid with constant density p. These assumptions are 
valid for the flow with high Reynolds and magnetic Reynolds numbers - 
such situation takes place in the system under consideration. The model in 
this case is simplifled considerably without removal any essential physical 
features. Also we note the paper P^ where the calculations of marginal 
stability for both viscous and inviscid fluids in the case of Couette flow are 
presented and they are shown to be almost identical. 

The fluid dynamics in our model is described by the equations of ideal 
incompressible MHD: 

^ + ( V-V) V = -VP + -1-{B-V)B, (2) 

V-y = 0, (3) 



— = Vx[VxB], (4) 

V-B = 0, (5) 

where P = {p + B • B/8'n')/p is the norinahzed total pressure. In equihbrium 
state {d/dt -^ 0) the fluid, placed in uniform magnetic field Bq = BqBz, 
rotates differentially with velocity Vq = rQ{r)e^f,. In general, the function 
Q{r) is arbitrary, since equilibrium can always be maintained by a radial 
pressure gradient, but we will consider here only the form of Q{r) defined in 
equation (P). 

For the purpose of linear stability analysis, it is advantageous to use a 
Lagrangian representation. We linearize the equations 0-® near the equi- 
librium state and introduce the displacement vector ^, which is considered 
to be a small quantity. The perturbations of velocity and magnetic field are 
then expressed in linear approximation as follows ^IGj : 

5V = ^ + (Fo-V)C-(^-V)Fo, (6) 

6B = Vxl^xBo]. (7) 

The linearized dynamics of ^ is obtained from equation (J2I): 

II + 2{Vo-V)^ + {Vo-Vn-{^-V){Vo-V)V,= (8) 

= -V6P+-^{6B-V)Bo + -^{Bo-V)6B. 
Airp Anp 

Besides, the equation (jH)) gives: 

V-| = 0. (9) 

The normal mode solutions to the system (jH))-© can be sought in the 
form: 

^(r,0,2,t) = ^(^)e^(™'^+^-— *), (10) 

SP{r,(j),z,t) = (5P(^)e*(™'^+'=^-^*), (11) 

where m (integer) and k are the azimuthal and axial wave-numbers, re- 
spectively, and u is the eigen-frequency. Since normal modes (jTUj). (fTT|) are 
generally complex, to find the perturbations of physical quantities one has 
to take either real or imaginary parts of righthand sides in equations (jHl), ^ 
and (HD). 



Substituting the normal mode representation (fTUj) into equations (jHl) and 
^ we find the perturbations: 



. . on 



5V ■■ 
5B = ikBoC 
while the system ©-© can be rewritten as 

or 



d{5P) 
dr 



{u\ — uj'^)iip — 2iVLuj^r 

{uj\ - uo^)^, = -tk6P, 
1 a(re 
r dr 



ira 



5P. 



r 



0, 



(12) 
(13) 

(14) 

(15) 
(16) 
(17) 



where we have introduced the shifted eigen-frequency, 

u) = uo — mQ, 
and the Alfven frequency, 

kBo 



Ua = kVA 



V^tFp 



The analogous derivation for more general case including toroidal magnetic 
field is given by jTTj. 

For computational convenience we introduce dimensionless quantities, 
taking as a unit of length the radius of the inner cylindrical boundary ri 
and as a unit of frequency the angular velocity at this radius Qi. Then we 
have: 

f = — , fi = 1, f2 = — , h = — , k = kri, 
ri ri ri 

Va ^ uj ^ uja r~ ^i~. 1 



Va 



where r2 is the radius of the outer cylindrical boundary. In the following 
consideration we will omit tildes to simplify notations. 

The system (IT ^ - (|T7jl can be reduced to one differential equation which 
in non-dimensional terms is: 



(cj^ — uj\)xu' 
rn^ + k'^x 



u 



u 



UJ, 



+ 



rau 



VI' - 



k^Q 



+ 



Ax ' m^ + k'^x \ m? + k'^x 

eQ{xQ' + Q) 



kW^Q^ 



+ 



(m^ + k'^x){uj'^ — C(j^) rn^ + k'^x 



+ 
0, (1^ 



where u = rC,r, x = r"^ and the prime denotes the derivative with respect to 
X. The further reduction of equation (fTHj) is possible if the function Q{r) of 
the form (^ is used. In this case Q oc 1/r^ oc 1/x and the last term in the 
equation is zero. 

Since we consider the channel with rigid walls the velocity components 
which are perpendicular to the boundary should vanish at that boundary. It 
gives us the following boundary conditions 

u{xi) = u{x2) = 0, (19) 

where Xi and X2 correspond to inner ri and outer r2 cylindrical boundaries 
respectively. Besides, we have to quantize axial wave-number k appropriately, 
namely, 

k = -^, (20) 

where Uz is integer and h is the height of the channel. 

The equation (|TH|l with boundary conditions ^T^ constitutes an eigen- 
value problem. The primary objective of our stability study is to determine 
the spectrum of eigen-frequencies u which has to be found as a part of the so- 
lution. Analyzing this eigenvalue problem one can conclude that its solution 
(and eigen-frequency spectrum as well) depends on five parameters: three 
wave-numbers - azimuthal m, axial k and radial rir (it will be introduced 
later), geometrical factor X2 and dimensionless Alfven velocity va- In this 
work we will focus on study of the dependencies of eigen-frequencies u on 
wave- numbers m, k and n^,. 

To solve the equation (fTHj) subject to the boundary conditions ffT^ we 
used a standard shooting method adapted for eigenvalue problem. For fixed 
values of problem parameters this method allows one to find the eigen- 
functions u{x) and corresponding eigen-frequencies u. We have performed 
calculations for the parameters of the experimental device proposed in Kur- 
chatov Institute (see Table Q). In numerical calculations the convenient nor- 
malization of eigen-function u'{xi) = 1 has been chosen. The following sec- 
tions represent the results of our calculations. 



3 Axisymmetric Perturbations, m = 

In the axisymmetric case with fl{x) = 1/x, equation (fT^ is reduced to 



u" + k\ 



_±+ ^i 



4a; ' x^{uj'^ — uj\y 



0. (21) 



Table 1: Parameters of experimental device. 



Parameter 


Value 


Dimensionless Value 


Inner radius, ri 


3 cm 


1 


Outer radius, r2 


15 cm 


5 


Height, h 


6 cm 


2 


Velocity at inner radius, vi 


32 m s-i 


1 


Alfven velocity, va 


8 ms"^ 


0.25 



Multiplying this equation by complex conjugate u* and integrating over the 
interval [xi,X2] we obtain the following variational form 



k^uol 



^2 

_ I /|2 _ _| |2 I 

Ax x^{uj'^ — u!\y 



u\ 



dx = 0. 



(22) 



XI 



Some conclusions can be made from its consideration: 



1. The expression (u;^ — uj 



2 \2 



a is a real positive number. Therefore, 



UJ 



UJ 



± a is real and eigen-frequencies uj are either pure real or 



pure imaginary. 

2. Eigen-frequency spectrum is symmetric with respect to both axis in 
the complex cu-plane. 

3. Eigen-f unctions u can always be chosen to be real. 

4. There is a degeneracy: each eigen-function u corresponds to four eigen- 



frequencies UJ. Two of them are always real, uj - 
means stability. Two others can be imaginary, uj 



±a/ci;^ + a^, which 
±a/w^ — a^, which 



means instability. This instability is known as the axisymmetric mag- 
net orot at ional instability (MRI of axisymmetric modes). 

Besides the azimuthal and axial wave-numbers, m and k, to describe 
eigen-frequencies completely we have to introduce the radial wave-number n^ 
which is defined in the case m = as a number of zeros of corresponding real 
eigen-function. The eigen-spectrum can be represented now by the formal 
disperse relation u = uj{m,k,nr) which relates each set of wave-numbers 
(m, k, Hr) with values of eigen-frequencies u. 

For m = and lowest axial wave-number k = 7r/2 (n^ = 1) the spec- 
trum of eigen-frequencies is shown in figure |21 As one can see, there is a 
pure imaginary eigen-frequency in the upper half-plane which corresponds 
to magneto-rotational instability. Also the spectrum has two accumulation 



points at uj = icu^i (Alfven resonances). With approaching these points the 
radial wave-number increases, i. e. the corresponding eigen-functions be- 
come oscillating. In this limit the WKB approximation to the equation (J2T| 
works well (see, for example, |12j). 

Eigen-spectrum shown in figure El is formed from coupled Alfven and slow 
magnetosonic modes. Unlike ^H] we did not obtain the fast magnetosonic 
modes since our consideration is restricted by incompressible MHD model. 

It should be noted here that the Alfven resonance lines. Re a; = ±uja, 
divide the complex tu-plane into 3 regions: 1) Re u; < —ua, 2) —uja <Re uj < 
lua, and 3) Re a; > uja- The "outer" modes belonging to regions 1 and 3 are 
always stable (the corresponding eigen- frequencies are real). The unstable 
modes can emerge only in the central region 2 as a result of coupling of 
originally stable "inner" modes. As we will show this tendency is valid for 
non-axisymmetric modes as well. 

In figure IHl the samples of axisymmetric eigen-functions with different 
radial wave-numbers rij. are shown. No conclusion about localization of these 
modes can be made in general. 

4 Non- Axisymmetric Perturbations, m ^ 

In the case of non-axisymmetric perturbations one has to solve the full equa- 
tion ()18|) in order to obtain the eigen-frequencies and make a judgment on 
the stability of the system. Some important properties of the eigen-frequency 
spectrum can be established prior to solving the problem: 

1. The eigen-frequencies of non-axisymmetric modes may assume gener- 
ally complex values. 

2. In the complex cu-plane the eigen-frequency spectrum is symmetric with 
respect to the line Im a; = 0, that is if u is an eigen-frequency then its 
complex conjugate u* is also eigen-frequency. The presence of complex 
conjugate pair in the eigenvalue spectrum indicates the instability of 
the equilibrium under consideration since one of these eigen-frequencies 
corresponds to the perturbation growing in time. 

3. Eigen-functions u are generally complex. 

4. There is no degeneracy in this case: each eigen-frequency corresponds 
to a different eigen-function. 

Unlike the axisymmetric case, in the case of m 7^ the Alfven resonances 
are not the lines any more. Instead they form two "resonance" zones in the 
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complex tu-plane: 

min {— uj A + 'mQ{x)) < lie uj < max {—ujyi + mQ{x)), (23) 

x£[xi, X2] xG[xi, X2] 

min (a;^ + mr2(x)) < Re a; < max {uA + mQ{x)), (24) 

x£[xi, xi\ x^\x\, X2] 

where the resonance condition Re cj = ±u;^ + mQ{x) is satisfied for some 
point X G [xi, X2]. As follows from our numerical analysis for each set 
of wave-numbers (m, k, rir) there are four different eigen-frequencies. Two 
of them are always real (stable) and situated just outside the "resonance" 
zones, one in the left outer region with 

Re u; < min {—oja + mVL{x)), 

x&[x-i, X2] 

another one in the right outer region with 

Re cu > max {oj a + 'mVL{x)) . 

x^[x-i, X2] 

The other two eigen-frequencies lie either in the "inner" region between the 
resonances 

max {—ujA + rnQ{x))<Reuj< min {uA + fnQ{x)), 

x£[xi, X2] x£[xi, X2] 

or in the "resonance" zones if these zones are overlapped. In former case the 
eigen-frequencies can be real or complex conjugate depending on the values 
of problem parameters. In latter case the eigen-frequencies are necessarily 
complex conjugate, that means the presence of instability. 

As a typical example of non-axisymmetric eigen-frequency spectrum the 
spectrum for m = 1 modes is shown in the figure IH The resonance zones 
in this case are overlapped and all the "inner" eigen-frequencies are complex 
conjugate (due to the resolution of the figure IH this fact is not so obvious 
for modes with the small imaginary parts of eigen-frequencies) . In figure El 
samples of these unstable eigen-modes are presented. 

Note that in the case of unstable mode with very small imaginary part 
of eigen-frequency a narrow singular layer is formed around the resonance 
point Xr defined by the condition 

Re u; = ±.uja + mVL{xr)- 

The eigen-function has the steep gradients in the neighborhood of this point 
(see Fig. |33,|3t). 
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For the unstable non-axisymmetric mode the definition of its radial wave- 
number Ur is not so obvious as for axisymmetric one because the non- 
axisymmetric eigen-function is complex. We resolve this problem as follows: 
changing m continuously from its initial value to m = we find the corre- 
sponding axisymmetric mode with known Ur, this rir is taken to be the radial 
wave-number of original non-axisymmetric mode. Although such definition 
is mathematically correct and non-ambiguous sometimes it leads to confusing 
results. For example, comparing the figures Eb andEb one can see that both 
real and imaginary parts of eigen-function with rir = 2 have less number of 
zeros than the corresponding parts of eigen-function with 71^ = 1. 

Our main numerical results are related to the study of MRI growth rate 
(or increment, 7 = Im w) as a function of wave-numbers m, k, riy. The figure 
ini demonstrate the dependence of instability growth rate on the azimuthal 
wave-number m {m is used here as continuous parameter) for eigen-modes 
with different radial wave-numbers rir. As one can see beginning with n^ = 1 
the system under consideration has a threshold of instability, rricr, besides 
the value of rricr is larger for larger n^. li m > rricr the eigen-mode with 
corresponding radial wave-number is always unstable, although the increment 
of instability decreases with growth of m. This fact contradicts to the local 
approach in which the unstable modes become stabilized with increasing 
values of m (see, for example, lUj). Also it is worth noting that the maximum 
possible value of the increment decreases with increase of radial wave-number 
rir- In this sense the most unstable mode is the mode with the lowest radial 
wave- number, rir = 0. 

The figure U\ demonstrates the dependence of the increment 7 on the 
azimuthal m and axial k wave-numbers for modes with n^ = 0. According 
to this dependence, for every value of m there exists a threshold value of 
axial wave- number kcr{m) after which the eigen-mode with n^ = is always 
stable. This threshold value increases with growth of m and, as we believe, 
approaches the asymptote kcr oc m (the explanation is in the next section). 

Dependencies of the increment on k for some m are shown in more details 
in figure |H1 In fact, this figure represents the sections of the figure [7| at 
corresponding m. 

Since the minimal axial wave-number k is determined by the height of the 
channel (see eq. Hd^ we conclude that choosing the height to be small enough 
(that is the minimal k is bigger than kcr) one can suppress the instabilities 
with low m. In particular, the situation is possible when there is instability 
for m = 1 mode but there is no instability for axisymmetric m = mode. 
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5 Marginal Stability 

As we have mentioned above if the resonance zones in the complex u -plane 
are overlapped all the inner modes become unstable. Although we checked 
this statement only for the small range of m (from to 3) we believe that 
this is also true for larger m. Our limitations in m are related to the computer 
accuracy: in fact the growth rates for large m are so small that they can not 
be resolved in the code even by double float representation. 

Accepting this statement as a hypothesis we obtain from conditions (|^H|l 
and ()24|) that all inner modes (with any radial wave- number n^) become 
unstable if 

max {mQ{x) — uja) = min {niQ{x) +^7^). 

xg[xi, X2] xe[xi, X2] 

Taking into account that xi = 1 and the angular velocity profile is Q{x) = 

1/x, we obtain 

m 

m — ujA = h oJAi 

X2 

which gives a critical value of axial wave-number A; as a function of m: 

fc,,(m) = ^fl-iy (25) 

It means that all "inner" modes (with any n^) become unstable if the axial 
wave-number k is less than the critical value defined by ()25|) . 

In reality the instability can occur for larger values of k, the actual 
marginal stability curves depend on the radial wave-number Ur- In the figure 
inithe marginal stability curves calculated for n^ = and Ur = 1 are shown. 
For comparison, the critical line kcr{m) is also plotted. It is quite clear that 
this line is asymptote for both presented curves in the limit m — »• cxd. Besides 
the critical line almost coincides with the marginal stability in the case of 
large radial wave-numbers, n,. ^ 1. 

Note that the system under consideration has the only unstable axisym- 
metric mode with radial wave- number n^ = 0. All axisymmetric modes with 
Ur > 1 are stable regardless of the value of k. From the figure El and expres- 
sion (^3)) it also follows that in this system there are always non-axisymmetric 
eigen-modes which are spectrally unstable due to magnetorotational insta- 
bility. 



6 Conclusions 

We have examined the spectral stability of electrically driven liquid metal 
flow in the model relevant to the new MRI experiment and presented a de- 
ll 



tailed numerical analysis of eigen-functions and eigen-frequencies inherent in 
this model. The obtained results show that the rotation profile given by the 
expression ((1} is always unstable due to MRI. For the chosen parameters of 
the experimental device the highest growth rate of instability 7 corresponds 
to the axisymmetric mode m = with the lowest radial wave- number Ur = 0. 
Consequently, in the experiment with electrically driven flow one should ex- 
pect the rapid excitation of axisymmetric magnetorotational instability. The 
modes with other wave-numbers are also presented but they are not so robust 
as the axisymmetric one. 

From our analysis it turns out that the growth rate of axisymmetric MRI 
is quite sensitive to the device parameters. For example, choosing the height 
of the channel to be small enough (to make the smallest possible k in the 
system bigger than its threshold value for m = 0) one can suppress the 
axisymmetric MRI while the unstable modes with m ^ will still exist. 

The model which we have considered here is, of course, idealized. We did 
not take into account the dissipative effects (viscosity and resistivity) as well 
as the deviations of the equilibrium velocity from the profile given by the 
equation (^. These effects can change the conditions for MRI excitation in 
the experiment; they will be a subject of our future study. 

We would like to acknowledge E. P. Velikhov for statement of this problem 
and V. P. Lakhin for useful discussions. This work is supported in part by 
NSERC Canada, NATO Collaboration Linkage Grant and Global Partners 
Award (University of Saskatchewan, Canada). 
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Figure 1: Geometry of the electrically driven flow for MRI study. 
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Figure 3: Eigen-functions in case m = 0, k = 7r/2. Corresponding radial 
wave-numbers and eigen-frequencies are (a) Ur = 0, u = ±0.2099i, ±0.5937; 
(b) nr = l,uj = ±0.1908, ± 0.5216; (c) rir = 2, lu = ±0.2625, ± 0.4894. 
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Figure 4: Eigen- frequency spectrum in case m = 1, k = 7r/2. Dashed lines 
represent real and imaginary axes, dotted lines correspond to boundaries of 
Alfven resonance zones. Radial wave-numbers n^ are shown for some unstable 
modes. 
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Figure 5: Eigen- functions in case m = 1, k = tt/2. Real (solid) and imaginary 
(dashed) parts are shown. Corresponding radial wave-numbers and eigen- 
frequencies are (a) Ur = 0, u = 0.5623 + 0.0791i; (b) n^. = 1, u; = 0.4173 + 
0.0073i; (c) nr = 2,uj = 0.3994 + 0.0093i 
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Figure 6: Growth rate dependence on azimuthal wave-number m in case 
k = 7i/2: (a) rir = 0; (b) rir = 1; (c) rir = 2. 
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Figure 7: Growth rate dependence on azimuthal m and axial k wave-numbers 
in case n^ = 0. Solid contours show the levels of function 7(m, k). The dashed 
line is the asymptote for the contour '~f{m, k) = 0. 




Figure 8: Growth rate dependence on axial wave-number k for different m in 
case fir = 0. Solid line is for m = 0, dashed - for 171 = 1, dotted - for m = 2. 
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Figure 9: Marginal stability in (/c, 'm)-plane for different radial wave-numbers 
n^ = 0. Solid line is for n^ = 0, dotted - for n^ = 1; dashed line is their 
asymptote. 
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